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Abstract 

A computation of the NRQCD (3 function both in the Lorentz gauge family and 
in the Coulomb gauge to one loop order endorses a velocity power counting scheme 
for dimensionally regularised NRQCD. In addition to the ultrasoft scale represented 
by bremsstrahlung gluons and the potential scale with Coulomb gluons and on-shell 
quarks, a soft regime is identified in which energies and momenta are of order Mv, 
gluons are on shell and the quark propagator becomes static. The instantaneous 
gluon propagator has a non-zero vacuum polarisation only because of contributions 
from this regime, irrespective of the gauge chosen. Rules are derived which allow 
one to read up from a given graph whether it is zero because of the homogeneity 
of dimensional regularisation. They also apply to threshold expansion and are used 
to prove that ultrasoft quarks with energy and momentum of order Mv 2 decouple 
from the theory. 
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1 Introduction 



Non-Relativistic QCD [I], [| takes advantage of the fact that when a heavy quark is nearly 
on shell and its energy is dominated by its mass M, the resulting non-relativistic system 
exhibits two small expansion parameters: the coupling constant g and the particle velocity 
v. The Coulomb interaction rules the level spacing in Charmonium and Bottomium 
because a s is small enough for perturbative calculations and the relative velocity of the 
quarks is v ~ a s (Mv) by virtue of the virial theorem, where the scale of g is set by the 
inverse Bohr radius Mv of the system. Albeit a s increases with decreasing Q 2 ~ (Mv ) 2 , 
a window between the relativistic perturbative and the confinement regime remains in 
which both a s and v is small (e.g. for Bottomium, a s (M^,v) ~ 0.35). Therefore, wave 
functions and potentials obtained by re-summation of ladder diagrams involving relevant 
couplings as v — > may be used to account for bound state physics, and calculations 
of production cross sections, hyperfine splittings, lifetimes, threshold properties etc. are 
much facilitated. The Abelian analogue, NRQED, to which all considerations in this 
article apply equally well, has also simplified precision calculations in Positronium. An 
incomplete list of recent references may include P|-||12||. 

The NRQCD Lagrangean in terms of the heavy quark (anti-quark) bi-spinors Q (Q) 
and gluons (D^ = <9 M + igA^) reads 

^NRQCD = Q t (i9o-^ )Q + ^Q ] D 2 Q + ^Q ] D A Q + ... 

+ H <$* ■ BQ + . . . - (CQ^aQ) ■ (Q ] aCQ) + ... (1.1) 

- J V'" + F£ v D 2 F fll/,a + ...+ £ GFix , 

where the coefficients q, di, e*, . . . encode the ultraviolet physics: All excitations with 
four-momenta of the order of M and higher are integrated out and give rise e.g. to four- 
point interactions between quarks {d\ ^ 0). The perturbative part of the coefficients 
can be determined by matching NRQCD matrix elements to their QCD counterparts 
in the regime where both theories are perturbative in the coupling constant. This has 
been performed to 0(M~ 3 ) At tree level, a Foldy-Wouthuysen transformation gives 
Ci = l,ei = 1, and loop corrections are down by powers of g, the most famous example 
being the coefficient for the Fermi term related to the anomalous magnetic moment of the 
electron, c 3 = 1 + ^ + . . .. Further coefficients to enter at one loop level are di = = 1. 
Lorentz invariance demands C\ = c-i = 1 to all orders. 

The Lagrangean ( |1.1[ ) consists of infinitely many terms constrained only by the symme- 
tries of the theory and is non-renormalisable. Predictive power is nonetheless established 
when only a finite number of terms contribute to a given order in the two expansion 
parameters^, g and v. Besides the heavy quark mass M, the typical binding energy and 
momentum scales in NRQCD are the non-relativistic kinetic energy Mv 2 and the mo- 
mentum Mv of the quark, which in Quarkonia appear as the energy and inverse size of 

-"^For clarity, the two will be distinguished in the following, although v ~ a s , as noted above. 
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the bound state (T|, |^|. Since for the smallest of the three scales a s (Mv 2 ) ft 1 in Bot- 
tomium and Charmonium, an expansion in g is justified only for the interactions taking 
place on scales Mv and higher in the real world, but one can imagine a world in which 
Mv 2 3> Aqcd- Toponium fulfils this requirement but decays mainly weakly and very fast 
so that QCD effects do not dominate. 

Because the effective Lagrangean does not exhibit the non-relativistic expansion pa- 
rameter v explicitly, a power counting scheme has to be established which determines 
uniquely which terms in the Lagrangean must be taken into account to render consistent 
calculations and predictive power to a given accuracy in v. It is at this point that NRQCD 
can serve as a "toy model" for nuclear physics Jl4| (although this grossly understates its 
valued): It will establish what the relevant kinematic regimes and infrared variables are 
in a theory with three (or more) separate scales, and it will demonstrate how to count 
powers of the non-relativistic expansion parameter, v. Recently, velocity power count- 
ing rules were established for a toy model \TE\ following Beneke and Smirnov's threshold 



expansion Jl6[, which has become very useful for understanding non-relativistic effective 
theories. The relevant energy and momentum regimes and low energy degrees of freedom 
were found for dimensionally regularised non-relativistic theories. This article presents 
the extension to NRQCD and as an application the calculation of the NRQCD (3 function. 
It will also be argued that the diagrammatic approach advocated here sheds a new light 
on threshold expansion |1| and helps simplify calculations. 

Both from a conceptual and a technical point of view the calculation of the NRQCD 
(3 function proves interesting. NRQCD is a well-defined field theory of quarks and gluons. 
One can therefore investigate its renormalisation group equations under the assumption 
that perturbation theory is applicable at all scales. For example, the running and mixing 
of the dimension six operators has been investigated to one loop order by Bauer and 
Manohar |~7|]. On the other hand, as NRQCD is the low energy limit of QCD, it must 
reproduce the QCD (3 function with Np light quarks below the scale M, but this has not 



been demonstrated so far. The problem with the power counting developed until now [14 



|T8| , was that gluons mediating the Coulomb interaction ("potential gluons") seemed 
to have a vanishing vacuum polarisation and hence a (3 function which was neither gauge 
invariant nor agreed with the QCD result, while the coupling of gluons on a much lower 
scale which mediate bremsstrahlung processes ran as expected. The apparent impasse 
which results is resolved in this article. 

Straightforward as the computation may thus seem with its outcome already antic- 
ipated, the prime goal of this article is not the result for the NRQCD (3 function; the 
objective is rather to shed more light on power counting in NRQCD and to provide a 
comparison between calculations in the effective and the full theory. It will show that 
the Lorentz gauge family (and indeed any standard gauge) is a legitimate gauge choice 
in NRQCD, not only the Coulomb gauge. It will prove that the recently discussed soft 
regime |T5| , [16], [2(| with quarks and gluons of four-momenta of order Mv is indispensable 
for NRQCD to describe the correct infrared limit of QCD. It will also test the validity 



2 The most striking difference to nuclear physics is that NRQCD has no exceptionally large scattering 
length to accommodate. 
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of the power counting proposed. In the Coulomb gauge, it will demonstrate that the 
renormalisation of the quark propagator and of the leading order quark gluon vertex are 
trivial, so that for the (3 function in this gauge, only a computation of the gluon vacuum 
polarisation is necessary. From a technical point of view, dimensional regularisation for 
non-covariant loop integrals proves to be a simple regulator choice. It allows to develop a 
set of simple diagrammatic rules to determine which graphs yield non-zero contributions 
in perturbation theory. This reduces the computational effort considerably as compared 
to other regularisation schemes in which power divergences are present like e.g. when us- 
ing a lattice cut-off. Because of the association of threshold expansion and NRQCD, such 
rules can be applied to the former, too. Finally, the calculation of the NRQCD (3 function 
presented here is - especially in the Feynman gauge - simpler than its QCD counterpart 
and shows some pedagogically intriguing aspects. 

The article is organised as follows: In Sect. ||], the velocity power counting of NRQCD 
is presented. The relevant regimes of NRQCD are identified (Sect. |2.1| ), extending the 
formalism of Luke and Savage PH by the soft regime of Beneke and Smirnov jl6|. The 



section then proposes the rescaling rules necessary for a Lagrangean with manifest veloc- 
ity power counting (Sect. |27|) and gives the vertex (Sect. |2^3|) and loop velocity power 
counting rules (Sect. ^4|) . Gauge invariance of the power counting is addressed in Sect. 
2~5| . The calculation of the NRQCD (3 function in the Lorentz gauges in Sect. starts by 



highlighting the intimate relation between threshold expansion and the proposed NRQCD 
power counting, which will additionally be summarised in the Conclusions. Diagrammatic 
rules developed in Sect. |3]2| help to facilitate the calculation considerably and are used 
in Sect. |3.3| to prove the absence of ultrasoft heavy quarks. As a tilly, the (3 function 
is computed in the Coulomb gauge in Sect. [3]7]. Summary and outlook conclude the 
article, together with an appendix containing exemplary calculations in non-covariant 
dimensional regularisation. 



2 Velocity Power Counting 



2.1 Regimes of NRQCD 



The first power counting rules for NRQCD were derived by Lepage et al. |H[ in the 
Coulomb gauge using the consistency of the equations of motion and a momentum cut- 
off. Only Coulomb gluons with typical momenta of order Mv were considered, taking 
into account retardation, but neglecting bremsstrahlung effects. The power counting also 
held to all orders only after all power divergences were subtracted. 

Simpler counting rules were proposed by Luke and Manohar for Coulomb interac- 
]14| |, and by Grinstein and Rothstein for bremsstrahlung processes |19 |. Luke and 



tions 

Savage |18|] united these schemes using dimensional regularisation, so that the tree level 
power counting is automatically preserved to all orders in perturbation theory up to log- 
arithmic corrections. They also extended the formalism to include the Lorentz gauges. 
Labelle proposed a similar scheme in time ordered perturbation theory. Beneke and 
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Smirnov [ HJ observed that the collinear divergence of the two gluon exchange contribution 
to Coulomb scattering between non-relativistic particles near threshold is not reproduced 
in this version of NRQCD because the soft low energy regime is not taken into account. 
A recent article [15j has shown that an extension of the work in |T3|, ISJ, |I^] using dimen- 
sional regularisation, as motivated from threshold expansion |yj , resolves this conflict in 
a toy model. The complete power counting scheme of NRQCD is presented briefly in the 
following (see also f20fl). In Sect. [371], the intimate relation between NRQCD as developed 
in recent years fl4|, [15], [Tj| and threshold expansion is addressed in more detail. 



The NRQCD propagators are read up from the Lagrangean Ql.l ) as 

i Num , ,, i Num 



Q 



T 



A" 



k 2 + ie 



(2.1) 



where T = p — M = ^ + ... is the kinetic energy of the quark. "Num" are numerators 
containing the appropriate colour, Dirac and flavour indices, all of which are unimportant 
for the considerations in this section and will be suppressed throughout this article. Gauge 
dependence of the gluon propagator for the power counting is addressed in Sect. [2~5| . 

Cuts and poles in scattering amplitudes close to threshold stem from saddle points of 
the loop integrand, corresponding to bound states and on-shell propagation of particles 
in intermediate states. They give rise to infrared divergences and in general dominate 
contributions to scattering amplitudes. With the two low energy scales at hand, and 
energies and momenta being of either scale, three regimes are identified in which either 
the quark or the gluon in ( |2.1| ) is on shell, as noted by Beneke and Smirnov 



soft regime: 
potential regime: 
ultrasoft regime: 



k 
T > 

h 



Mv 2 



Mv , 
\p\ ~ Mv 



(2.2) 



\k\ ~ Mv 2 



Note that the scale M has been integrated out when deriving the non-relativistic La- 
grangean (|1.1|) , so that a "hard" regime with on shell gluons and quarks of four-momentum 
of order M cannot be considered. 

What is the particle content of each of the three regimes? Ultrasoft gluons A 1 ^ are 
emitted as bremsstrahlung or from excited states in the bound system, and hence are 
physical. Soft gluons Ag do not describe bremsstrahlung: Because in- and outgoing 
quarks Q p are close to their mass shell, they have an energy of order Mv 2 . Therefore, 
overall energy conservation forbids all processes with outgoing soft gluons but without 
ingoing ones, and vice versa, as their energy is of order Mv. Finally, gluons which change 
the quark momenta but keep them close to their mass shell relate the (instantaneous) 
Coulomb interaction: 

Ap 1 : k ~ Mv 2 , |jfe| ~ Mv (2.3) 

When a soft gluon A^ couples to a potential quark Q p , the outgoing quark is far off 
its mass shell and carries energy and momentum of order Mv. Therefore, consistency 
requires the existence of quarks in the soft regime as well, 



Q s : T ~ \p\ ~ Mv 



(2.4) 
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With these five fields Q S} Q p} Ag, Ag, representing quarks and gluons in the three differ- 
ent non-relativistic regimes, soft, potential and ultrasoft, NRQCD becomes self-consistent. 
These fields are the infrared-relevant degrees of freedom representing one and the same 
non-relativistic particle in the respective kinematic regimes and came naturally by identi- 
fying all possible particle poles in the non-relativistic propagators. Their interactions are 
fixed by the non-relativistic Lagrangean 

( PP . Section p| will prove that a hypothetical 
ultrasoft quark, created e.g. by the radiation of a potential gluon off a potential quark, 
decouples completely from the theory. On the other hand, Sect. ||] will demonstrate that 
in order to obtain the correct result for the NRQCD (3 function, all fields listed above 
in the three regimes have to be accounted for. Therefore, the particle content presented 
above is not only consistent but both minimal and complete. 

In order to guarantee that there is no overlap between interactions and particles in 
different regimes, the regular isat ion scheme must finally be chosen such that the expan- 
sion around one saddle point in the loop integral as performed above does not obtain 
any contribution from other regimes, represented by other saddle points in the loop in- 
tegrals. One might use an energy and momentum cut-off separating the soft from the 
potential, and the potential from the ultrasoft regime, but the integrals encountered can 
in general not be performed analytically and in closed form. Furthermore, by introducing 
another, artificial scale, cut-off regularisation jeopardises power counting and symmetries 
(gauge, chiral, . . . ). Power divergences occur when the (un-physical) cut-off is removed 
in intermediate steps but not in the final, physical result, complicating the computation. 
In contradistinction, using dimensional regularisation after the saddle point expansion 
preserves power counting because its homogeneity guarantees that contributions from dif- 
ferent saddle points and regimes do not overlap. (A simple example is given in Ref. ||15|| .) 
Homogeneity will also be essential when developing diagrammatic rules to classify graphs 
as zero in Sect. |3.2| and when showing the decoupling of the ultrasoft quark in Sect. |3.3|. 



2.2 Rescaling Rules and Propagators 



In order to establish explicit velocity power counting in the NRQCD Lagrangean, one 
rescales the space-time coordinates such that typical momenta in either regime are di- 
mensionless, as proposed by Luke and Manohar [14| for the potential regime and by 
Grinstein and Rothstein fOl for the ultrasoft one: 



soft 
potential 
ultrasoft 



x 
x 

X 



{Mv) 
(Mv) 



(2.5) 



(Mv)~ l T s , 

{Mv 2 )- l T u , x = (Mv 2 )' 1 X u . 

For the propagator terms in the NRQCD Lagrangean to be normalised as order v°, one 
sets for the representatives of the gluons in the three regimes |15|, 

soft: A»(x,t) = (Mv) A£(X S ,T B ) , 
potential: A£{x, t) = (Mv?) A%(X S , T u ) , (2.6) 
ultrasoft: A£{x,t) = {Mv 2 ) A£(X n ,T n ) , 
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and for the quark representatives 

soft: Q B (x,t) = {Mv)i Q S (X S ,T S ) , (2.7) 
potential: Q p (x,t) = (Mo)? Q P (X S ,T U ) . 

The rescaled free quark Lagrangean reads then 

soft: d 3 X s dT s Qj(id + V -d 2 )Q s , (2.8) 

potential: d 3 X s dT u Qj(id + - <9 2 ) Q p . (2.9) 

Here, as in the following, the positions of the fields have been left out whenever they 
coincide with the rescaled variables of the volume element. Derivatives are to be taken 
with respect to the rescaled variables of the volume element unless otherwise stated. In 
order to maintain velocity power counting, corrections of order v or higher must be treated 
as insertions, so that one reads up the (un-rescaled) quark propagators and insertion as 

(T,p) i ( T 'P) £2 

soft: Q s : » : , )( : -i -f— = P.C.h 1 ) ,(2.10) 

yes T + ie 2M 

(T,p) j 

potential: Q p : » : ^ — . (2-H) 

The soft quark becomes static because T ~ Mti ^> ^ ~ Mu 2 allows one to expand the 
non-relativistic propagator ( |2.1| ) and to treat the momentum term as an insertion. 

Throughout this article, the symbol P.C(t> n ) denotes at which order in the velocity 
power counting a certain term or diagram contributes. In the approach presented, each 
line in a loop diagram counts as P.C(t>°) because the strength of the propagator has been 
set to unity in the rescaled Lagrangean ( |2.8| / [2.9| ). The rescaled Lagrangean measures the 
strength of the insertion relative to the propagator. Therefore, the soft quark insertion 
(|2.10|) does not count as which scales like v 2 , but as %m/T = P.C.(f). In contradis- 



tinction to the approach in [16], one will therefore not obtain the absolute order in v at 
which a given graph contributes, but the relative order in v between graphs or vertices is 
read off more easily and asserted correctly. 

The gauge fixing term was included in the NRQCD Lagrangean ( |1.1| ) since the decompo- 
sition of the Lagrangean into a free and an interacting part is gauge dependent. Because 
of the difference between canonical and physical momentum, it is important to specify 
the gauge before identifying to which order in v a certain regime in the Lagrangean con- 
tributes, as will be seen shortly. In the following, the Feynman rules for the Lorentz 
gauges and for the Coulomb gauge are derived explicitly. Sect. will comment on the 
gauge invariance of the procedure. 

The rescaled free gluon Lagrangean in the Lorentz gauges reads 



1 



2 



1 



soft: d s X s dT s - A% <9V " (1 " -W» A" , (2-12) 



a 
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potential: d 3 X s dT u - Alt g^(v 2 d 2 - & 



(2.13) 



(1 - -)OV)<9o + S^d^vS^do + 5 vi di)\A u u , 



ultrasoft: d 3 X u dT u \ A* [<9 V - (1 - -)9^J ^ , 



while in the Coulomb gauge 



soft: 



1 



(B 2 - dl)5 l3 - d % d 3 



A 



■3,8 ' 



d d X s dT s - A itS 

l3v JT> ^ f A 5^ A i A i&X ;);) .,2 ;) 2. 



potential: crX s dT u - [A, P <9 A) iP + A, P (9 <% - d t dj - v d Q 5ij)A 3 ^ 
ultrasoft: d 3 X s dT s ^ A, u [(<9 2 - d 2 )5 i3 - ^d,]^ . 



(2.14) 

(2.15) 
(2.16) 
(2.17) 



The (un-rescaled) Coulomb and Lorentz gauge propagators are therefore given as {SH(k) 



soft: A{! : WWWWWV: 



Coulomb gauge Lorentz gauges 
i^tr(^) iGpv(k) 



potential: A PtQ : 
A p : 

ultrasoft: At : 



k,A o 



fc,A 



P + ie 
— i 

-jfe 2 + ie 

i jg(j) 
-k 2 + ie 

i g(g 
k 2 + ie 



k 2 + ie ' 
— i 

-fc 2 + ie ' 
i + (1 - a>-^ 



(2.18) 



fc 2 +ie 



-/c 2 + ie 
A; 2 + ie 



The potential gluon becomes instantaneous in both gauges as expected for the particle 
mediating the Coulomb interaction. Insertions are necessary only in the potential regime: 



Coulomb gauge 



kAo o 
----X--- : 

1 k > A 3 

----X--- : 

i 

----X--- : 



none 



i fcg 5ij 



P.CV 



none 



Lorentz gauges 

-iM = P . C > 2 ) 

a 

ik 2 5 %J = Y.C.{v 2 ) 



-i (1 - -) = P-C» 

a 



(2.19) 



The Lorentz gauge propagators and insertions for the potential and ultrasoft regimes were 
first given in ||18| . Especially for the Feynman gauge a = 1 {Q^ v {k) = g^ u ), Lorentz and 
Coulomb gauges in the potential regime differ only by insertions, i.e. at higher order in v. 
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As seen from fl2.15| - [2T7l ), the choice of the Coulomb gauge makes A instantaneous 
to all orders in v, and hence it contributes in the potential regime, only. Since in this 
gauge, A solely mediates the instantaneous Coulomb potential (physical fields are trans- 
verse by virtue of GauB' law), this result was to be expected. The field A p is associated 
with retardation effects like spin-orbit coupling and the Darwin term in In the 

Coulomb gauge, the advantages of A Q only contributing in the potential regime and of its 
propagator having no insertions or admixtures with the vector components of the gauge 
field, are balanced by the fact that a more lengthy and cumbersome renormalisation seems 
necessary. The calculation of the NRQCD f3 function will prove that the Lorentz gauges 
are a legitimate gauge choice in NRQCD, although the number of diagrams appears to 
be larger because of a larger number of vertices, as will be seen now. 



2.3 Vertex Rules 



By experience, particles in the various regimes couple: On-shell (potential) quarks radiate 
bremsstrahlung (ultrasoft) gluons. In general, one must allow all couplings between the 
various regimes which obey "scale conservation": Both energies and momenta must be 
conserved within each regime to the order in v one works. This will exclude for example 
the coupling of two potential quarks (T ~ Mv 2 ) to one soft gluon (go ~ Mv), but not 
to two soft gluons via the Q^A ■ AQ term of the Lagrangean (IT). For a more rigorous 
derivation of this rule see Sect. |Q| . 

As an example, consider a bremsstrahlung-like process: the radiation of a soft scalar 
gluon T s off a soft quark, resulting in a potential quark with four momentum (T^p 1 ). 
The rescaled interaction Lagrangean with its Hermitean conjugate reads 



d 3 A s dT s 



-gv° Qj(X B ,T s )A B , {X s ,T s )Q p (X s ,vT s ) + H.c. 



(2.20) 



Note that the scaling regime of the volume element is set by the particle with the highest 
momentum and energy. Therefore, one obtains the power counting of a given vertex in 
a specific regime. As the external quarks are on shell, i.e. in the potential regime, the 
power counting rules of subgraphs whose scaling rules have been determined in the soft 
regime have to be transferred to the potential regime. This step is postponed to the next 
section. The strength of the example vertex above is read up as v in the soft regime. 
This is clearly also the strength of the Hermitean conjugate vertex. 

The interaction Lagrangean is non-local in the rescaled variables, and in order to 
maintain velocity power counting, Q p (X s ,vT s ) has to be expanded about Q p (A s ,0) in 
powers of v. The Feynman rule for this vertex is after inverting the rescaling 




-ig (2tt) 4 5 {3) (p -p' + q)x 
exp(-T^ q^)) S(T- 



¥.C.{e v ) . 



(2.21) 



The Hermitean conjugate vertex is built analogously. As the potential quark has a much 
smaller energy than either of the soft particles, it can - by the uncertainty relation - 



9 



not resolve the precise time at which the soft quark emits or absorbs the soft gluon. 
This "temporal" multipole expansion comes technically from the different scaling of x 
and t in the three regimes. The multipole expansion symbolised in the power counting 
by e v corresponds term by term to an expansion in v and should be truncated at the 
desired order. In general, the coupling between particles of different regimes will not 
be point-like but contain multipole expansions for the particle belonging to the weaker 
kinematic regime. That the coupling of potential quarks to ultrasoft gluons requires a 
momentum multipole expansion as in atomic physics has been observed by Grinstein and 
Rothstein pf, and by Labelle J22 . 

Amongst the fields introduced, six scale conserving interactions are allowed within and 
between the various regimes for any coupling of one gluon to two quarks. Because only in 
the Coulomb gauge, A is a field in the potential regime only, its propagator always being 
instantaneous, only the first two interactions exist for the scalar coupling —gQ^A Q of 
table [I] in this gauge choice. The v counting for the lowest order quark-gluon interactions 
from this vertex is presented in table |l|. 



Table 1: Velocity power counting and vertices for the interaction Lagrangean —gQ^A^Q. 
In the Coulomb gauge, only the first two diagrams exist. The last line indicates the field 
for which an energy or momentum multipole expansion has to be performed. 



Vertex 














I 

I A 
1 


I 
I 
I 










P.C. 


1 


i 

V 2 


v° 


v° 


v° 


V 


multipole 


none 


A»(x) 


A° u (x) 


none 


Q P (t) 





Note that - although both describing interactions with physical gluons - soft and 
ultrasoft couplings occur at different orders in v. On the level of the vertex rules, an 
overlap of different regimes resulting in double counting is prevented by the fact that in 
addition to most of the propagators, all vertices are distinct because of different multipole 
expansion rules. 

This approach may be compared to power counting in loop diagrams as proposed in 
Beneke and Smirnov's threshold expansion |16| where the strength of all scalar gluon 
interactions is v° because the scalar vertex Q^AqQ of table [I] does not contain derivatives. 
Here, one can easily see that the Coulomb interaction is the only relevant coupling as 
v — ► and that it scales like t> _1//2 . This follows immediately from the rescaling rules 
proposed. Also, the suppression of bremsstrahlung processes relative to the Coulomb 
interaction is evident from table |I| In threshold expansion, these features are established 
by considering scattering processes with an arbitrary number of loops. 

Velocity power counting for other vertices is again obtained by rescaling and multipole 
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expansion. For later reference, the counting rules for the coupling of one and two gluons, 
either minimally or via the Fermi term, and the rules for the three gluon interaction are 
displayed in tables ||, |3] and f|. In the second, couplings between on shell gluons in the 
same regime are all of order v° = c = 1 as expected for relativistic particles. 

In the minimal coupling term — 1| Q^d ■ AQ, the derivative acts on both the gluon 
and the quark field. Because the quark is either soft or potential, its derivative from 
QU • dQ is rescaled as d -> {Mv) d s . The same holds when the derivative acts on a soft 
or potential gluon. But both for the one gluon part of the Fermi interaction Q^a ■ BQ 
and for the term Q\d ■ A)Q of the minimal coupling, the derivative acts on an ultrasoft 
gluon field and must hence be rescaled as d — > {Mv 2 ) d u . That this part of the vertex 
coupling is one power in v weaker than the part where the derivative acts on the quark, 
is recorded in parentheses in table The v rules for the coupling of one gluon to the 
quark via the Fermi term are identical to those sub-dominant contributions to the minimal 
coupling vertex for ultrasoft gluons, and to the dominant contributions in the other cases. 
A similar provision is made in table [| but is not necessary in table |j. 

Table 2: Velocity power counting and vertices for the interaction — 1| Q^d ■ AQ. When 
the spatial derivative acts on an ultrasoft gluon field, the power counting in parentheses 
is the sub-dominant contribution as explained in the text. It coincides with the power 
counting for the one gluon component of the Fermi vertex. 

















Vertex 


i 

1 A 








1 
1 
1 




P.C. 


l 

V 2 


v (v 2 ) 


V 


V 


3 

V 2 
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Table 3: Velocity power counting and vertices for the gluonic interaction Lagrangean 
2f abc (dnA^ — d u A^) A b ^ A c ' u . Sub-dominant contributions in parentheses. 
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2.4 Loop Rules 

As hinted upon above, the velocity power counting is not yet complete. The propagators 
were constructed such that they count as order v in each regime in a Feynman diagram. 
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Table 4: Velocity power counting and vertices for the coupling of two gluons to a quark 
from the minimal coupling or the Fermi term in the NRQCD Lagrangean fll.lp . 
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As one sees from the volume element used in ( 2.20| ), the vertex rules for the soft regime 
count powers of v with respect to the soft regime. There, one hence retrieves the velocity 
power counting of Heavy Quark Effective Theory |23|, |24| (HQET) , in which the interac- 
tions between one heavy (and hence static) and one or several light quarks are described. 
HQET counts inverse powers of mass in the Lagrangean, but because in the soft regime 
Mv ~ const., the two approaches are actually equivalent. Now, HQET becomes a sub-set 
of NRQCD, complemented by interactions between soft (HQET) and potential or ultra- 
soft particles. This has already been used by Manohar Jl3| to facilitate and unify the 
matching of the NRQCD/HQET Lagrangean to QCD. 

In NRQCD with on-shell quarks as initial and final states, the soft regime can occur 
only inside loops. Let us define as a "soft blob" each connected graph containing (soft) 
loops which is obtained when all potential and ultrasoft lines are cut. If the diagram 
contains even one soft particle, scale conservation ensures that there is at least one loop 
which consists of only soft particles, be they quarks or gluons, and that it is part of the 
soft blob. Inside the soft blob, the power counting for the vertices is performed in the 
soft regime and has therefore to be transfered to the potential regime. Since soft loop 
momenta scale like [d 4 /c s ] ~ t> 4 while potential ones like [d 4 fc p ] ~ v 5 , each soft blob is 
enhanced by an additional factor -. 

Consider for example the graphs of Fig. [3]: Using the Lorentz gauges, vertex power 
counting gives that the leading contribution is from the exchange of two potential gluons, 
coupled via Q^AqQ. There are four such vertices, so the diagram is 0(g A ) P.C.(i; -2 ) from 
table |1|. The next two diagrams are O(g G )P.C.(v ) and 0(g 2 ) P.C.(t> x ) from the vertex 
power counting in the soft regime, but another factor - must be included because there 
is one soft blob in each diagram. This way, the v counting of the soft regime is moved 
to the potential one. The intermediate couplings in the second diagram take place in the 
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Figure 1: Power counting with soft loops. The loops in the second and third diagram 
obtain an inverse power of v, the last diagram of v 2 in addition to the power counting 
following from the vertex rules. 



soft regime and hence are counted in that regime. After cutting all potential and ultrasoft 
lines in the last diagram, two soft blobs are separated by the propagation of two potential 
quarks. The graph is 0(g 14 ) P.C(v °) from the vertices, and the loop counting gives a 
factor — j. Each soft blob contributes at least four powers of g, but only one inverse power 
of v ~ g 2 . Power counting is preserved. These velocity power counting rules in loops are 
also verified in explicit calculations of exemplary graphs, [TE\ and App. |A.2| . 



There is no similar rule for ultrasoft loops: In the absence of ultrasoft quarks (see Sect. 
|3.3j ), the internal ultrasoft gluon couples ultimately to a particle in the potential or soft 
regime. Those vertices are automatically counted in this stronger regime, while couplings 
between ultrasoft particles are counted in the weakest regime. No "ultrasoft blobs" can 
therefore be isolated by cutting all potential and soft lines. It is hence ultimately scale 
conservation which forbids a non-trivial loop counting rule for the ultrasoft regime. 

Irrespective of the gauge chosen, there is only one relevant quark-gluon coupling at 
tree level in the renormalisation group approach (i.e. only one which dominates at zero 
velocity): As expected, it is the Q P A p>0 Q p coupling providing the binding (table [TJ). The 
potential gluon ladders must therefore be re-summed to all orders to yield the - Coulomb 
potential. Diagrams higher order in v are corrections. In the Coulomb gauge, all other 
couplings and insertions are irrelevant, while in general, there are three marginal couplings: 
Q P A n fiQ p , QlA^oQs and QlA St0 Q p . Because of the additional factor v^ 1 per soft blob, 
graphs containing the latter two couplings can indeed be relevant and contribute as v — > 
(e.g. the second graph in Fig. |]). Eventually, retardation effects from A 1 ^ will become 
weaker than contributions from the soft regime. 

One finally turns to the inclusion of other relativistic particles. In the same way as 
NRQCD replaces the physical gluon with one representative per regime, any light (rel- 
ativistic) particle has to be tripled. There are therefore three ghost fields n s ,n p ,r] u with 
the same rescaling rules as the gluon fields ( |2.(j| ). Fermions require more thought: In 
the real world, the kinetic energy of the b quark in Bottomium is compatible to the 
strange quark mass, M\>v 2 ~ m s . For the sake of simplicity, this article assumes all light 
particles to have masses very much smaller than any other scale, m q <C Mv 2 , so that 
the relativistic particles can be treated as massless to lowest order and the denominators 
of the light particle propagators are identical to the ones of the gluon in the respec- 
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tive regimes. The Dirac spinors representing light quarks scale in the three regimes as 
ip s ~ (Mf)a ~ ipp,ip n ~ (Mt) 2 )2, i.e. similar to the heavy quark but including its ul- 
trasoft counterpart (see also Sect. |3.3|) . The number of vertices per term in the rescaled 
interaction Lagrangean increases as shown above. Quark-ghost couplings and heavy-light 
quark couplings will appear in the Lagrangean (|1 . 1|) at 0(g 4 ) to restore unitarity and can 
be obtained by matching NRQCD to QCD. 

With rescaling, multipole expansion and loop counting, the velocity power counting rules 
are established. To witness, the rescaling rules (|2.5| / |2.6| / |2.7|) provide an efficient and 



well-defined way to arrive at an NRQCD power counting: After rescaling, one reads up 
the order in v at which any term in the Lagrangean contributes in each of the three 
kinematic regimes and performs the appropriate multipole expansions. This establishes 
the Feynman rules for NRQCD and HQET simultaneously, and classifies the strength of 
the vertices in the regime of the particle with highest energy and momentum. Finally, the 
rescaling is inverted, introducing one un-rescaled gluon and quark field as the infrared- 
relevant degrees of freedom for each kinematic regime. To obtain the power counting 
for a certain graph, loop counting is taken into account to transfer the strength of a 
soft blob into the potential regime. Computations are performed most naturally in the 
original, dimensionful variables. The rescaling rules are only needed to establish the power 
counting, but in order to maintain it, one is of course not allowed to re-sum the multipole 
expansions in the un-rescaled variables. It is also interesting to note that there is no 
choice but to assign one and the same coupling strength g to each interaction. Different 
couplings for one vertex in different regimes are not allowed. This is to be expected, as 
the fields in the various regimes are representatives of one and the same non-relativistic 
particle, whose interactions are fixed by the non-relativistic Lagrangean (|1 . 1| ) . 



2.5 Gauge Invar iance 

The Coulomb gauge d ■ A = is a natural choice in NRQED because in it, static charges 
do not radiate. In NRQCD, gluons carry colour and hence the Coulomb gauge does 



not have this advantage. Luke and Savage [18| showed how to establish explicit velocity 



power counting in the Lorentz gauges of NRQCD, too. The set of gauges applicable to 
NRQCD can be extended further: The classification of the three kinematic regimes fl2.2|) 
itself relied only on the typical excitation energy and momentum, and hence on gauge 
invariant quantities, and on the form of the denominator in the propagators fl2.1|) which 
is unchanged in any order of perturbation theory^. The perturbative quark propagator 
is gauge independent. Gauge fixing will introduce gauge dependent denominators mul- 
tiplying the gauge independent denominators in the perturbative gluon propagators. As 
an example that this will in general not change the identification of the soft and ultra- 
soft regimes from poles in the gluon propagators, consider the generalised axial gauge 



3 



Non-perturbatively, the propagators are not of the form (^1) because of confinement and the absence 



of coloured states, so that the power counting presented here may break down. 
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(n 2 = —1, a arbitrary), in which 



k 2 + ie 



k u n u + njz v 1 + ak 2 

9r> - 



(2.22) 



k ■ n (n ■ ky 

The additional denominators n-k introduce no new combinations of the two infrared scales 
Mv and Mv 2 for which the gluon propagator has a pole. Therefore, the decomposition 
into the three regimes ( |2.2| ) remains unchanged, as do the rescaling properties of the fields 
and interactions, (gTg/gg/g77|) and tables |T| to f§ 



The standard gauges (axial, Weyl, Lorentz, Coulomb) will therefore all show the same 
power counting and vertex rules quoted above. Details of the gluon propagator and 
its insertions look different in different regimes and gauges, and some gauges will not 
exhibit certain vertices, insertions and representatives, e.g. the Coulomb gauge is unique 
in having A contribute only in the potential regime. Only when the gauge dependent 
denominator introduces a new regime has the power counting to be modified. This is for 
example the case in a gauge with denominator k$ — M 2 k 2 , in which an exceptional regime 
(k ~ Mv , \k\ ~ Mv 2 ) enters. Rescaling is then performed including the new regime. 

As is well known, the Weyl gauge A = wants a constraint quantisation: GauB' 
law D ■ E = gQ^Q, generating local gauge transformations, is the equation of motion 
derived by varying the Lagrangean with respect to Aq. When Aq = 0, GauB' law will 
not be recovered as an equation of motion. In order to restore local gauge invariance, 
a projector onto states obeying GauB' law has therefore to be inserted into the path 
integral. Resolving GauB' law explicitly, the Lagrangean of the Coulomb interaction has 
the structure 

^coulomb = / d 3 y g 2 Q^Q(x) Q{x - y) Q ] Q(y) , (2.23) 

where Q{x — y) is an appropriate instantaneous Green's function to GauB' law with dimen- 
sion [mass] 1 , e.g. in QED T+fl. Using the rescaling rules (p|/|2^/|2^) , one derives the 



Coulomb interaction between two quarks to be P.C.(t> x ) as in the other gauges. Without 
constraint quantisation, the longitudinal component of the vector gauge field mediating 



the Coulomb interaction in Weyl gauge QED, A 1 , has a static propagator 



3 The NRQCD (3 Function 

This chapter presents the computation of the perturbative part of the NRQCD (3 function 
to order g 3 , v° in the MS scheme using the vertex and rescaling rules derived in the 
previous section. Initially, the Lorentz gauges are not only chosen because they allow a 
less cumbersome renormalisation than e.g. the Coulomb gauge, but also to demonstrate 
that in the final result, the gauge parameter a drops out. The Slavnov-Taylor identities 
will be shown to be fulfilled, and the Lorentz gauges are a legitimate gauge choice in 
NRQCD. As a by-product, Sect. will calculate the NRQCD (3 function to lowest order 
in the Coulomb gauge. 

By construction, NRQCD and QCD must agree in the infrared limit, and especially 
in the structure of collinear (infrared) divergences. Matching proved that both the soft 
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quark and the soft gluon are indispensable to reproduce the correct structure of collinear 
divergences in a toy model given by Beneke and Smirnov [TJJ and confirmed the proposed 



counting rules [HJ . The calculation of the NRQCD (3 function will manifest the relevance 



of the soft regime beyond infrared divergences, and it will endorse the power counting 
further. Deriving the lowest order (3 function in the following, the relation between this 



power counting and threshold expansion |16[ will be exemplified, and rules will be devel- 



oped which allow one to determine from the structure of a diagram whether it is zero or 
not to all orders in v. 

NRQCD with Np light quarks must reproduce the (3 function of QCD below the scale M, 



Pqcd — - 



3 3 



(4vr 



(3.1) 



to lowest order for the gauge group SU(iV) and renormalised coupling g^. This means 
especially that the renormalised coupling strengths of all interactions steming from ex- 
panding the same term in the Lagrangean in the various regimes are the same, except that 
they have to be taken at different scales: For the interaction Lagrangean —gQ^A Q, the 
six vertices of table [I] all couple with the same un-renormalised strength. Although the 
number of vertices is increased in the approach presented here, the number of independent 
couplings is not. Indeed, one major result of this chapter will be that renormalisation will 
not destroy this symmetry because the renormalisation constants agree in all regimes. 

One can obtain the (3 function from the ghost-gluon coupling, and the ghost and gluon 
self energy. Because the relativistic sector of NRQCD is identical to the one of QCD as 
seen in Sect. |2j| , the calculation proceeds in that case like in QCD and yields the same 
result. In this article, the (3 function is calculated by studying the heavy quark sector, 
namely the renormalised coupling Qr^o.rQr between the heavy quark and the scalar 
gluon. This term yields immediately the coupling renormalisation. All other vertices are 
suppressed by at least one power of v. 

Because the integrals to be performed are not Lorentz invariant, standard formulae for 



dimensional regularisation often do not apply. Therefore, Ref. [15] used split dimensional 
regularisation, which was introduced by Leibbrandt and Williams |25| to cure the problems 
arising from pinch singularities in non-covariant gauges. It treats the temporal and spatial 
components of the loop integrations on an equal footing by regularising the energy and 
momentum integration separately, J -0^ = J |^ , a -»• 1, d -»• 4 |26|, Chap. 

4.1]. Useful formulae for the computation of the (3 function are given in App. [A.l] . 



3.1 NRQCD and Threshold Expansion: Potential Quark Self- 
Energy 



Threshold expansion [[HJ and NRQCD make use of very similar basic techniques but 
different formulations, as this Section outlines. In order to clarify the relation, consider 
the quark self-energy diagram lowest order in g which couples the scalar gluon and the 
quark in NRQCD. To facilitate the presentation initially, the Feynman gauge a = 1 is 
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chosen in which scalar and vector gauge fields do not mix when propagating. Without 
the scaling rules for the various regimes, one obtains 

A ,q 




: iEp(r , ? ) = ( - i#CF /^-^_^_, (3.2) 

with Cpdij = (t a t a )ij = N l^ l 5ij the Casimir operator of the fundamental representation of 
SU(iV). Threshold expansion identifies the loop momentum of QCD to belong to a hard 
regime q ~ M or to either of the three regimes (|2.2| ) and expands the integrand about the 
various saddle points, i.e. about the values of the loop-momentum q where poles occur, like 



the NRQCD classification in Sect. |2.1| . The hard regime subsumes the relativistic effects 
contained in the coefficients q, d iy of the NRQCD Lagrangean ( |1.1| ) and hence provides 
the matching between NRQCD and QCD |16|]. In deriving the NRQCD Lagrangean ( |1 . 1|) , 



hard momenta have been integrated out, so that the ultraviolet behaviour of the above 
loop integral is arbitrary and it is not necessary to take the poles at q ~ M in ( j3.2|) into 
account. Indeed, to be consistent, they should rather be discarded. 

Let the incident quark be on shell, i.e. have a four-momentum (T p ,p) in the potential 
regime. The integral decomposes then into approximations about three saddle points: 
When q is soft, the gluon propagator has a pole, and one can expand the quark propagator 
in powers of T p /q ~ v and /l° ~ v - ^ ne <l uar k becomes static: 

1 ' X +-ifc-^±|)!U + ... (3.3) 



1o + T p - ^ go go V p 2M J g 



iE p (T,p) = (-is) 2 CV / — V ^ ] (3.4) 

^ '^soft ^ y > J (2-K) d q 2 Q - g 2 + ie q + ie ^ \ q + ie ,l ; 



The NRQCD power counting proceeds on the level of the Lagrangean instead of the 
Feynman diagrams with corresponding diagrams for loop momenta in the soft regime, 



+ * X > + ••• (3.5) 

\ 2 r f ddq ~ [ 1 (1 ■ (P + v) 2 1 i-r 1 \ 

1 l9) ° F J (2nYqi - q> + ie q + ie \ L * 2M q + ie + Up q + ie + ' ' ') ' 

and recovers order by order in v the result of threshold expansion. The intermediate 
soft quark is static, and the higher order terms in threshold expansion are interpreted 
as insertions into the soft quark propagator or as resulting from the energy multipole 
expansion at the Q S A£Q P vertex. In fact, using the equations of motion, a temporal 
multipole expansion may be re-written such that the energy becomes conserved at the 
vertex. Now, both soft and potential or ultrasoft energies are present in the propagators, 
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making it necessary to expand it in ultrasoft and potential energies. An example would 
be to restate the vertex ( |2.21| ) as 



: -ig(2n)U(T-T; + q )5^(p-p' + q) , (3.6) 
and the soft propagator to contain insertions P.C.(f) for potential energies Tp 

77 E M • (3-7) 



j — : — i — l_ 



9o + ie^ V go 

The same can be shown for the momentum-non-conserving vertices, too. 

When the loop momentum q is potential, threshold expansion picks the pole of the 
quark propagator, and the gluon propagator is expanded as if go -C |g|, corresponding to 
NRQCD diagrams with insertions in the gluon propagator: 

i 

— + ... (3.8) 

iS p( T '^Ln tial = { ~' l9)2 ° F I^Y ^TYe n ? (l) T + g -W + i e 

The correspondence is again complete for ultrasoft q, and the higher order terms of thresh- 
old expansion come from the spatial multipole expansion in NRQCD: 



iS p (T,p) = (3.9) 

PV ultrasoft V ' 



2M 



-■\*r f ddq ~ [ i ^ 

l9} ° F J (2nY a 2 n-Q 2 + ie T + a»-£, + \c 



(27r) d g 2 -g 2 + ie T + g - g + ie to V T + ?o - + « 



Threshold expansion and NRQCD concur indeed in all diagrams to be calculated for the 
(3 function. Further comparison of the two approaches is postponed to the Conclusions. 



3.2 Some Diagrammatic Rules 

In general, the number of possible diagrams in power counted NRQCD is considerably 
larger than in the original theory because there are at least six vertices per interaction 
allowed by scale conservation, cf. tables [3] to f| for interactions between three and four 
particles. In the preceding sub-section, three diagrams were drawn for the lowest order 
contribution to the potential quark self energy. The analogue problem in threshold ex- 
pansion is that each loop momentum can lie in either of the three regimes, so that the 
number of expansions for an n loop diagram grows like 3 n . 
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The important step to cut down the computational effort in either approach is to note 



that as an immediate consequence of the axioms of dimensional regularisation (e.g. |p6 
Chap. 4.1 and 4.2]), integrals without scales set by external momenta or energies in the 
denominator vanish because of homogeneity, 

d d a 

1 q a = . (3.10) 



(2tt) c 

For example, applying this theorem to the potential quark self energy diagrams, contri- 



butions from soft (3A) and potential ( |3.8|) loop momenta are zero to all orders in the 
velocity expansion. For the potential gluon contribution ( |3.8j ), this is seen by shifting the 
regularised loop integral T + q — > q before using (|3.10| ) for the energy integral. 



The importance of ( p.lOj ) has already been alluded to in threshold expansion [1~6|, and 
here a more thorough and formal treatment is presented. 



It is this theorem fl3.10| ) which renders most of the diagrams in NRQCD and threshold 



expansion calculations zero, and it is usually not even necessary to consider the whole 
diagram but only a sub- set of it. The following set of rules is helpful for reducing the 
numbers of diagrams to be dealt with in the calculation of the (3 function. 

In order to establish these rules, recall that all graphs vanish which contain a sub-graph 
zero by the rules developed. The routing of the loop four-momentum inside a diagram 
is arbitrary as always in dimensional regularisation so that having regularised the loop 
integration, all loop four-momenta can be shifted at will like in ordinary integration. 
Since they have identical denominators (Sect. |2.3| ), gluon lines in any rule can be replaced 
by any relativistic particle, e.g. ghosts and light quarks. Finally and most importantly, 
numerators are unimportant to determine whether an integral is scale-less following ( |3.10| ). 
As multipole expansion and insertions do not change the denominators of the propagators, 
a diagram is therefore zero to all orders when it is zero because of ( |3.10| ) at leading order 
in v . This result is also insensitive to the gauge chosen and to the specific vertex involved. 
For example, because the numerators play no role in rendering (O / |3.8| ) zero using ( 5TTD ), 



the same will hold for any graph with the same diagrammatic representation but different 
interactions, e.g. also for graphs in which one or both vertices are replaced with the 
minimal coupled vector fields Q^d ■ AQ or the Fermi interaction Q^a ■ BQ. 

As a first rule, consider a potential gluon between two soft particles, represented by the 
overlay of double and zigzag lines in Fig. [| (a). This sub-diagram is zero when the loop 
energy go is integrated over: Assign the potential loop momentum q to the instantaneous 
gluon with a propagator denominator — q 2 ( |2.18|) . As the coupling of soft particles to 
potential ones is energy non-conserving, the denominators of the soft particles coupling to 
do not contain q Q ~ Mv 2 , either. Therefore, qo does not occur in any denominator, and 
the dimensionally regularised loop integral over q is zero from Q3.10D . The rule is extended 



by noting (Fig. |2| (b)) that overall energy conservation allows for the potential gluon line 
to be attached with an arbitrary number of vertices into which an arbitrary number of 
potential or ultrasoft gluons or quarks with collective (potential) four-momentum I enters. 
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(T,p)T T(TV') 




(d) 

Figure 2: The first rule as proven step by step in the text, (a) Simple case; (b) a 
more complicated sub-diagram; (c,d) generalisations. The overlay of double and zigzag 
line stands for any number of arbitrary soft particles, the triple line for any number of 
potential or ultrasoft particles, all entering at the same vertex. The blob represents vertex 
and propagator dressings. No time ordering is implied in the way the diagrams are drawn. 

The denominators of the potential gluon between two soft particles still do not contain 
the loop energy, and the go integral vanishes again, Fig. |2| (c). With this four-momentum 
routing, I can be an external or loop momentum. 

One can finally dress the vertices and propagators without changing the argument 
thanks to scale conservation. Also, when several other potential or ultrasoft particles 
couple to the same soft particle in addition to the potential gluon, a temporal multipole 
expansion is necessary for each of them. Therefore, a potential gluon can never connect 
two soft particles, irrespective of the number of particles coupling to the same vertices as 
the potential gluon, and irrespective of other particles coupling to the potential gluon in 
the intermediate state, Fig. |2] (d): Diagrams with unbroken potential gluon lines between 
soft particles vanish. 

As predicted, the only vertex property which enters is that the coupling between 
particles of different regimes necessitates a multipole expansion which can be read up 
from the diagram directly. The proof therefore extends to massive, relativistic potential 
particles, like a light quark. It eliminates for example the diagram corresponding to a 
non-Abelian vertex correction with a denominator 

( T,p) (T',p >) 

* ZZz Zz/ q : Denom. = T' f ((T - T'f - (p - p' + qf) . (3.11) 
A similar rule can be proven for the potential quark between soft particles, see Fig. |3| (a) 
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for its bare and (b) for its dressed form. The vertices are again energy non- conserving, so 




(b) 



(c) 



Figure 3: The second rule in its bare (a) and dressed (b) version; (c) the generalisation 
analogous to Fig. $(c) fails. Conventions as in Fig. |^. 

->2 

that go enters only in the denominator of the potential quark in the combination go — 2&J- 

-*2 

Shifting Qo—2M^ %, the g integral again does not contain a scale in the denominator. In 
contradistinction to the previous case, this rule cannot be extended to include potential or 
ultrasoft momenta coupling to the blob, Fig. [3] (c). For example, a bremsstrahlung gluon 
with external, ultrasoft four-momentum / renders the potential quark denominators in 

->2 -"2 

the propagators as (go — 2m)(Qo + ~ 2m)- Now, l provides the necessary scale for the g 
integration. As an example of the rule, the Abelian vertex contribution with denominator 



(T,P) 



: Denom. 



(3.12) 



is zero, but the last of the graphs in Fig. [I] is not. 



As the coupling of an ultrasoft gluon to a soft particle conserves neither energy nor 
momentum, one may disconnect the ultrasoft gluon from such a vertex, Fig. ^ (a). If after 




Figure 4: (a) Soft-to-ultrasoft vertices are cut off in a third rule; (b) an example. Con- 
ventions as in Fig. 

cutting, a line with loop momenta to be integrated over becomes completely disconnected 
from the graph, the resulting tadpole makes the graph vanish, Fig. [| (b). If only one of 
the legs becomes dis-attached, the diagram can be non-zero, as the following diagram 
for the Abelian vertex correction to the Q s A s0 Q p coupling demonstrates: 



(T,P) 



(T',p') 



^ : Denom. = T q 2 \ T - q - 



2M 



(3.13) 
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See also Sect. 3J; and the computation of this graph in App. |A.2 . 

When several ultrasoft gluons couple to a soft particle at the same point, cutting does 
usually not result in a tadpole, either. For example, the gluon-gluon vertex obtained by 
cutting the two-gluon two-quark vertex 





^ : Denom. = (k - q) 2 q V 



(3.14) 



might be represented by a small blob as a reminder that it is neither energy nor momentum 
conserving. In the "tadpole" , a scale is set by the external ultrasoft four-momentum k. 

Loops made out of one soft quark and one soft gluon are restricted by another rule. In 





(b) 



Figure 5: The fourth rule in its (a) bare and (b) extended version, where the rule in Fig. 
[| (a) was used for attaching ultrasoft gluons to the soft gluon. Conventions as in Fig. |^. 

Fig. [| (a), the denominators are q 2 . — q 2 for the soft gluon and qo for the soft quark. 
Therefore, the diagram is without scale. The generalisation, Fig. [5] (b), eliminates all soft 
contributions to the potential quark self energy ( |3.5|) and to Abelian vertex corrections 
between potential and ultrasoft particles. Note that the light particle exchanged is mass- 
less or has a mass considerably smaller than Mv 2 . A relativistic particle mass of order 
Mv 2 or bigger provides a scale in the relativistic propagator, invalidating this rule. Also, 
as soon as a potential particle couples to the soft gluon, the diagram is non-zero. An 
example is the non- Abelian vertex correction to the Q^A£Q p coupling with denominator 

(T,P) (T',p f ) 




: Denom. = q q 2 (q 2 , — (p 



p' + qf 



(3.15) 



To summarise, the principle underlying the rules is simple: One identifies a potentially 
vanishing sub-set of the entire diagram and assigns the external and loop momenta to it, 
taking into account multipole expansions. Then, the denominators (i.e. inverse propaga- 
tors) are written out at lowest order in v . If by shifting integration variables, one arrives 
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at an energy or momentum integral without scale, the diagram vanishes to all orders. 
Along this line, further rules are easily established. 



Two more rules are more straightforwardly derived from the propagator properties of 
NRQCD directly. The first one cuts down the number of diagrams including potential 
gluons. For the heavy quarks, there is no distinction between Feynman and retarded 
propagators in NRQCD (|2.10|/pTTT ) because antiparticle propagation has been eliminated 



by the field transformation from the relativistic to the non-relativistic Lagrangean. The 
fields Q s and Q p propagate hence only time-like into the forward light cone. Still, Feyn- 
man's perturbation theory becomes more convenient than the time-ordered formalism, as 
less diagrams have to be calculated. Soft and ultrasoft gluons propagate light-like, and 
potential gluons instantaneously. Therefore, any diagram which - when time ordered - 
involves quarks propagating outside the forward light cone are zero, as are diagrams which 
cannot be drawn with both ends of a potential gluon being space-like separated. This 
"forwardness" excludes all non-zero contributions from the potential gluon contribution 
to the quark self energy in any regime, e.g. ( |3.8|) . 

Another rule is that there must at least be one pole on each side of the real go ax is 
for a diagram to be non-zero. This will e.g. not be the case when the loop momentum q 
only occurs in one quark propagator besides potential gluon propagators. Likewise, the 
contribution to the Abelian vertex correction 
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Denom. = \q, - — + iej \q, — + iej (T - [p - q) 

(3.16) 

has only poles below the real axis. Closing the contour of the go integration above it, the 
diagram is zero. This can also be proven using dimensional regularisation. 

In order to demonstrate that the rules established above are useful in reducing the number 
of graphs, Fig. ^ shows that only one of the four scale-conserving graphs contributing to 
the one loop soft quark self-energy survives their filter. The classification of the only 
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(a) (b) (c) (d) 

Figure 6: Out of the four scale-conserving one loop contributions to the soft quark self 
energy, only the first one, (a), survives, (b) is zero because it contains a sub-graph zero by 
Fig. (a); (c) vanishes because of Fig. ^ (c) or because of forwardness; (d) is eliminated 
by cutting off the ultrasoft gluons (Fig. [| (a)). The resulting tadpole is zero. 

non-zero diagrams for Abelian and non- Abelian one-loop corrections to one gluon vertices 
addressed later, Figs. |], || and |K| may serve as another example. 
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In conclusion, the homogeneity of dimensional regularisation ( |3.10 ) can be translated 
into diagrammatic rules to systematically identify graphs which vanish to all orders in 
threshold expansion or velocity expansion in NRQCD. The concept is sensitive only to 
the multipole expansion of the vertices involved, not to their precise nature. It is also 
independent of insertions and gauge chosen. Only the denominators of the graph lowest 
order in v have to be looked at. 



3.3 Theorems from Diagrammatic Rules 

Diagrammatic rules are quite effective to prove decoupling theorems to all orders in per- 
turbation theory. "Scale conservation" entered in Sect. [2.3| as a constraint on both the 
overall energy and momentum of a vertex because both have to be conserved order by 
order in v. Now, one may resort to a more formal argument: If one rescales the fields 
in the Lagrangean in all possible combinations, diagrams containing scale non-conserving 
vertices are eliminated because they correspond to temporal or spatial tadpoles. As an 
example, consider a soft particle coupling to two potential ones. The dashed-solid line 
representing any potential particle, the multipole expansion is 

(T,p) (T',f) 



':UAo oc (2tt) 4 6®(p -p' + q) [exp ((T - T') A) S(q )} . (3.17) 

The temporal tadpole of the go integration makes all diagrams vanish in which this vertex 
is embedded such that q is a loop momentum. If q Q is an external energy, its being zero 
contradicts the assumption that it is of order Mv. Therefore, diagrams containing scale 
non-conserving vertices are zero. In threshold expansion, this follows automatically from 
the observation that not the individual vertices but the loop momentum as a whole is 



expanded, see Sect. 3.1 



Another application is the decoupling of an ultrasoft heavy quark from NRQCD. When 



the particle content of NRQCD in the three different regimes was outlined in Sect. \L1 
the bremsstrahlung gluon was assumed to be the only ultrasoft field, and an ultrasoft 
quark Q u was not necessary to make the theory self-consistent. An ultrasoft quark might 
still be produced by the radiation of a potential gluon off a potential quark, but it is not 
unavoidable. In contradistinction, a soft quark is indispensable because a soft, on shell 
gluon excites a potential quark to a kinetic energy of order Mv: The potential quark 
coupling to a soft gluon becomes necessarily soft. Taking advantage of the homogeneity 



( 3.10 ) of dimensional regularisation, one can show that a hypothetical ultrasoft heavy 
quark decouples completely from the theory: All graphs containing it are zero. This fact 
was already remarked upon by Beneke and Smirnov [16| , and the formal proof using the 



NRQCD techniques proceeds as follows. 

Of all properties of the ultrasoft quark, the only ones needed are that its propagator 
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(denoted by a dotted line) is static because T ~ Mv 2 , \p\ ~ Mv 2 , 



Qu-- 
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T + ie 



(3.11 



and that it couples momentum non-conserving to all but ultrasoft particles. These fea- 
tures can of course be derived from the rescaling rules of an ultrasoft quark, Q n (x,t) = 



(Mv 2 ) 2 Q n (X u , T u ), in full analogy to Sect. |2~2] . If it exists, Q u can only enter in internal 
lines. Fermion number conservation dictates that it is produced and annihilated in a 
vertex into which at least one soft or potential quark enters. 

The simplest sub-diagram containing an ultrasoft quark is depicted in Fig. [7] (a) and 
vanishes because there is no loop momentum q in any denominator. This is still true 
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Figure 7: Decoupling of the heavy, ultrasoft quark, (a) primitive diagram; (b) extension 
to the case of ultrasoft particles coupling to Q u ; (c) general vertex involving potential and 
ultrasoft particles; (d) ultrasoft quarks decouple from any graph containing potential and 
ultrasoft particles. Conventions as in Fig. 0. 

when an arbitrary number of ultrasoft particles couple to Q u , Fig. |7| (b), and when such a 
coupling occurs repetitively. Now one considers couplings to potential particles like in Fig. 
(c). The momentum multipole expansion necessary for all ultrasoft particles separately 
in such a vertex disallows q to be transferred to a denominator. Therefore, any diagram 
containing an ultrasoft quark coupled to potential or ultrasoft particles is zero because of 
the spatial tadpole in the q integration, Fig. [7] (d). Because couplings of ultrasoft particles 
to soft ones involve the same momentum multipole expansion as couplings to potential 
ones, the proof extends to all graphs containing ultrasoft quarks: Every diagram with 
ultrasoft quarks is zero. Recursive application dresses all vertices and propagators. 
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3.4 The Quark Self-Energy 

The rules cover the lowest order contributions to the potential quark self energy as dis- 
cussed in Sect. |Q| . The graph in (|3.5|) is zero because of Fig. |!| the one in (|3T8| ) because 



of forwardness. The only remaining graph is (|3.9| ), whose lowest order contribution is 
P.C.(f°) by the power counting of table |l|. Indeed, corrections to particle propagators, 
quark or gluon, must count as v since the power counting was constructed such that free 
particle propagators in each regime scale like v ( |2.6|/f2.7| ) and renormalisation corrections 
to propagators must be of the same order. Using the integral (|A.4|) of App. [A.l , the pole 
in four dimensions is extracted as 



^ P (T,p)U = ^^-C F ^T[e] ( T "|^)+ P -C>) , (3-19) 



where e = 2 — | and /i is the renormalisation parameter. The potential quark propagator 
is recovered, confirming that different regimes do not mix under renormalisation. The 
non-relativistic quark propagator does not need renormalisation of the heavy quark mass. 
The quark wave function renormalisation is hence in the MS scheme 



, ^ P r[e] . (3.20) 
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For the soft quark self energy, the only diagram surviving the diagrammatic filter (Fig. 
10) is at lowest order in the power counting using the Feynman rules and 




i£ s (T,p)U = Cf ^ r[e] T . (3.21) 

That the soft quark renormalisation is the same as for the potential quark is not surprising 
as both stem from the same un-expanded NRQCD diagram ( |3.2j ). 

3.5 The Vacuum Polarisation 

In contradistinction to the self energy for Q p which does not receive contributions from 
soft particles, the vacuum polarisation of the potential gluon is non-zero only because of 
the presence of soft gluons, independent of the gauge chosen. The scale-conserving graphs 
with gluon loops are 




V- -( >- . (3.22) 



The integral in the second graph does not contain q in the denominator and hence 
is zero; the third diagram vanishes because it couples a potential gluon to a light-like 
particle twice. Soft gluons are therefore indispensable to provide the gluons mediating 
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the Coulomb interaction with a non-zero vacuum polarisation, and hence to have a running 
(3 function. As the gluons in the soft and ultrasoft regimes are on-shell particles and must 
run in NRQCD as in QCD, there is no reason to expect the potential gluons to freeze 
out in perturbation theory. The only non-zero ghost and light quark contributions come 
analogously from soft ghost and soft light quark propagation in the loop. 

The rescaling rules of table [| give the three gluon vertex to count as v%, and an 
additional v _1 from the loop power counting makes the graph count as v at leading 
order. Again, this is expected not only because gluons are relativistic particles, but also 
because one wants to renormalise a propagator which is P.C(v°) in the power counting. 



Without the loop rule of Sect. |2.4j , the graph would be of order v 1 , and the power counting 
would predict that its contribution would vanish as v — > so that the potential gluon 
propagator at rest would reduce to the bare one, a clearly unacceptable conclusion. The 



explicit computation ( |3.24 ) will also show that this graph is P.C.( 



The soft and ultrasoft gluon receive their vacuum polarisations from loops with on shell 
gluons, P.C.(v°): 



o 



(3.23) 



All other diagrams vanish by the rules developed above. Both contributions are identical 
to the ordinary QCD result, as no insertions or expansions enter. The gluon vacuum 
polarisation from gluon, ghost and light fermion contributions is therefore in the soft and 
ultrasoft regime [E7|, eq. (2.5.132)] 
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(3.24) 



Because the potential gluon vacuum polarisation does not contain insertions or multipole 
expansions in the internal lines, either, the QCD result can be taken and expanded in 
powers of the external energy fco "C \k\. As the infinite part of U(k 2 ) does not contain k, 
the only change will be that to lowest order in v, the part guaranteeing transversality of 
the gluon becomes 



{kfjty, - k 2 g^) -> (SpiSvjkikj + k 2 g^ v ) + P.C.(f) . 



(3.25) 



Renormalisation therefore keeps the contributions from the three regimes separate and 
the potential gluon propagator transversal up to higher order in v. For all three regimes, 
the gluon wave function renormalisation is in the end the one of QCD |27], eq. (2.5.135)], 
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3.6 The Vertex Correction and NRQCD (3 Function 

Since they probe only the relativistic sector of the theory, the renormalisation constants 
for gluons, ghosts and other light particles are the same in QCD and NRQCD. The quark 
wave function renormalisation is computed in the non-relativistic sector and for a non- 
relativistic bi-spinor rather than a relativistic Dirac spinor, so that it is not surprising 



that the result ( ^20l) differs from its QCD counterpart Z 2 QCD = 1 - 9 -^r C F a T[e] |7 



eq. (2.5.139)] even in the dependence on a, although the MS scheme was used in both 
cases. As both Lagrangeans are gauge invariant and agree in the light particle sector, 
the Slavnov-Taylor identities of QCD must also hold in NRQCD. Therefore, the renor- 
malisation Z\ of the quark-gluon vertex Q^AqQ can be inferred from Z 2 ,Z 3 and the 
three-gluon renormalisation Zi t9 of QCD (and NRQCD) not to be the same as in QCD 

where Z? CB = 1 - &f (C F a + N 2±a) T[e] 0, eq. (2.5.145)]: 



^ = ^ = 1 + |^(Cp(3-o:)- N 3 -±^) T[s] ^Z? CD . (3.27) 

In the following, one will identify as the cause that all non-zero contributions from both 
Abelian and non-Abelian vertex corrections to the scalar gluon vertex in NRQCD are 
different from their QCD values: Most notably, the non-Abelian vertex only provides 
gauge parameter corrections. As a by-product, the topologies of all one loop corrections 
to interactions of one gluon with a quark will be found and their leading order velocity 
power counting determined. 

Which vertices may be encountered in the graphs leading in vl For every vertex, the 
Q^AqQ coupling is at least one order stronger in v than any other coupling (tables [l] 
and 0). Any contribution is therefore suppressed in which scalar gluons in a vertex are 
replaced with vector gluons, coupled either minimally or via the Fermi term fff Q ^ ■ BQ 
in ( |1 . 1| ) . Moreover, the Fermi term couples one or two gluons to the quark spin. When 
only one Fermi interaction is found in any vertex correction, parity conservation requires 
the strongest possible correction involving A as outgoing particle to be of the form of 
the spin-orbit term, i.e. proportional to Q^a ■ (D x E)Q. This cannot correct the leading 
scalar gluon vertex. Two Fermi interactions are even more suppressed because vector 
couplings are weaker than scalar couplings. 

Each Abelian correction graph (Fig. starts off at the same order as the bare graph 
whose vertex correction it presents (table [I]) when only scalar gluons couple to the quark. 
The only exception is the last diagram in the first row of Fig. |[ which by the power 
counting is P.C.(v x ) while the bare QlA s0 Q p vertex is P.C.(v°) from table [I], as confirmed 
by an explicit computation of this graph in App. [A.2| . This diagram does therefore not 



enter in the computation of the vertex correction at lowest order. 

The non-Abelian corrections to the Q^A Q Q vertex at order g 3 fall into two categories: 
The topology of the first class of diagrams, Fig. |9], is analogous to the one of the non- 
Abelian vertex correction in QCD. In the Feynman gauge a = 1, its leading order con- 
tribution involves two scalar gluons and one vector gluon because in that gauge, different 
components of the gluon field do not mix when propagating, (|2.18|) , and the coupling of 
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Figure 8: The Abelian vertex corrections of order g 3 from the interactions Q^A Q which 
survive the diagrammatic filter. The leading order power counting is obtained when all 
vertices are the leading order scalar gluon interactions. After applying the rules of Sect. 



cL2\ , the seven diagrams drawn are the only survivors of 22 scale conserving graphs. 



three scalar gluons is forbidden at order g by the Lorentz structure. Again because for 
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Figure 9: The non- Abelian vertex corrections of order g 3 which survive the diagrammatic 
filter. The soft blob in the second diagram in the second row provides an additional factor 
- by loop power counting. The power counting indicated is the leading contribution in the 
Feynman gauge, when the three gluon vertex is the standard QCD one and both one vector 
and one scalar gluon couples minimally to the quark. In the Lorentz gauges, two scalar 
gluons couple to the quark, making all diagrams a factor - stronger and hence contribute 
at the same order as the Abelian corrections in Fig. [|. Out of 25 scale-conserving diagrams, 
only eight survive the diagrammatic filter. 

every vertex, the Q^AqQ coupling is stronger in v than any other coupling, the NRQCD 
analogue to the non- Abelian QCD vertex correction to the scalar gluon vertex is sub- 
leading in the Feynman gauge. The resulting power counting is reported in Fig. In the 
generalised Lorentz gauges, the scalar and vector components of the gauge field mix in 
the gluon propagation, with the amount of mixing proportional to (1 — a) ( |2.18| ). A scalar 
gluon emitted from a quark will therefore partially turn into a vector gluon, which then 
enters a three gluon vertex. The non-Abelian vertex corrections now have a non-trivial 
Lorentz structure and evaluate to be proportional to (1 — a) and (1 — a) 2 , as is confirmed 
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in App. |A.2| , ( A.15| ). A factor - is gained with respect to the power counting with one 
vector gluon coupling as given in Fig. ^, making the diagrams of the same, leading order 
as the Abelian vertex corrections. The non-Abelian graphs can consequently be seen as 
merely providing gauge corrections to the Abelian vertex result. The last diagram in the 
first row and the next-to-last in the second row of Fig. |9| are sub-leading and vanish using 
the equations of motion at lowest order for the outgoing potential quark, see App. |A.2 . 



There is no comparable graph in QCD to the second category of non-Abelian vertex 
corrections in NRQCD at order g 3 . Figure lists the topologies of all non-zero corrections 
to one gluon vertices involving one gluon loop. Again, one convinces oneself easily from 
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Figure 10: The only non-zero diagrams of the vertex corrections of order g 3 from the Fermi 



term proportional to C3 in the NRQCD Lagrangean ( |Q} ). The soft blob in the next-to-last 
diagram provides an additional factor - by loop power counting. The leading order power 
counting is obtained when the three gluon vertex is the standard QCD one and the gluons 
couple to the quark via the Fermi term. Out of 24 scale-conserving diagrams, only six 
survive the diagrammatic filter. 



tables |3] and |] that these diagrams are at least one order in v weaker than the Abelian 
vertex corrections, even before the precise vertex structure is specified. In NRQED, this 
type of diagrams occurs first at order g 5 from the e<i term in the Lagrangean ( |1.1| ) because 
a three photon coupling does not occur earlier. In QCD, the lowest order three gluon 
coupling is totally antisymmetric in colour space, so that its contraction with the totally 
symmetric expression in colour space from the vertex j^Q^A 2 Q vanishes. The Fermi 
interaction has a colour antisymmetric two gluon part because Bf = ^e^ k (djAl — d^A a - + 
gf ahc A b j A c k ) ) but is not a possible source of an order g 3 vertex correction in Fig. [10] to the 
minimal coupling of scalar gluons and quarks because the spin structures do not match. 

Hence, the non-Abelian nature of the gauge field enters into the NRQCD (3 function 
only via the gluon vacuum polarisation and as what may be called gauge correction 
represented by non-Abelian vertex corrections. Each of the six quark gluon couplings 
between the various regimes is renormalised by one and only one non-zero Abelian and 
non-Abelian vertex correction. 

Since the classification of diagrams following the rules of Sect. |3]2| is insensitive to the 
precise nature of the vertex coupling, the list leads as a by-product to a classification of 
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all non-zero one loop graphs in which two quarks and one gluon couple. The same is true 
for the quark self energy and vacuum polarisation graphs. The v power counting of the 
diagrams in Figs. |], |9| and is specific to the interactions chosen, but the fact that all 
other graphs with the same topology are zero is not. 

All Abelian graphs lead to the same vertex vertex correction, as demonstrated in App. 
A.2j at exemplary graphs: 



- i g R ^ 1* r Abelian | pole = (C F - *L) t« T[e] (3.28) 

Likewise for the non- Abelian vertex corrections: 
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coincides therefore with the one expected from the Slavnov-Taylor identities ( 3.271) . 



The NRQCD (3 function is computed from the scale invariance of the bare coupling 

S = 2 R/ /-^= (3.31) 
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To the order calculated, one may "match" the relativistic and non-relativistic (3 functions 
to infer that indeed to lowest order in g and v, 

QCD NRQCD / Q Q o\ 

9r =9r ■ (3-33) 

To match NRQCD to QCD via the renormalisation group behaviour of its couplings and 
operators is of course considerably more complicated than to match matrix elements. 
Nonetheless, it demonstrates that matching relations survive renormalisation. The renor- 
malised NRQCD coupling is as expected the same in all regimes, gauge invariant, and 
runs simultaneously in all regimes. The NRQCD (3 function is independent of the gauge 
parameter a. The Lorentz gauges are a legitimate gauge choice in NRQCD, obeying the 
Slavnov-Taylor identities. Wave function and vertex renormalisations differ in QCD and 
NRQCD, while the /3 functions do not. 

In dimensionally regularised NRQCD, the vacuum polarisation receives its sole con- 
tribution from the propagation of on-shell relativistic particles, i.e. of physical soft or 
ultrasoft gluons and soft or ultrasoft ghosts and light quarks. Potential gluons, i.e. such 
mediating the Coulomb interaction, do not give rise to a gluon renormalisation, nor do 
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they contribute to any other renormalisation at lowest order in v, (|3.19|) and Figs. |6|, ||] 
and || This observation applies to all standard gauges, including physical ones like the 
Coulomb gauge. It is traced back to the homogeneity of dimensional regularisation (13.10Q , 



which in turn is well known to be connected to a cancellation of ultraviolet and infrared 
divergences in scale-less integrals, e.g. [^7], p. 172 f.]. In a cut-off regularisation, the situ- 
ation is different as massless tadpoles require both infrared and ultraviolet regularisation 
resulting in logarithms, invalidating (|3.10|) and all diagrammatic rules, e.g. |28], Chap. 



18.5]. In such a regularisation, the Coulomb gluon vacuum polarisation can therefore 
indeed receive its main contribution from Coulomb gluons in the intermediate state. 

To close this section, mind that this derivation of the NRQCD (5 function is slightly 
simpler than its QCD counterpart but may be applied with equal generality. Except in the 
determination of the gluon vacuum polarisation, the Lorentz structure of the vertices and 
propagators did not enter. The Dirac algebra is not needed except when light fermions 
are included, and then the manipulations are straightforward. The vertex corrections in 
the Feynman gauge reduce to the computation of one graph: the Abelian vertex which 
involves neither Lorentz nor Dirac indices. Non-Abelian vertex corrections are easily 
extracted, see App. |A.2j . The presentation demonstrates also intuitively that quarks with 
mass M higher than the scale at which the j3 function is computed freeze out and are not 
to be included in the light quark number Np of ( |3.32|) . In a world with Nf quarks, the 
QCD (3 function can hence be found from NRQCD with Nf = Np "light" quarks and one 
fictitious quark with a mass much bigger than any real quark mass. 



3.7 The Coulomb Gauge (3 Function 

With the gauge independent classification of all one-loop self-energy, vacuum polarisation 
and vertex correction diagrams completed, it is straightforward to calculate the (3 function 
in the Coulomb gauge from the renormalisation of the scalar gluon field. In order to 
classify the contributing diagrams, one combines the notion that in the Coulomb gauge, 
A is only potential (Sect. |2.2| ) with the knowledge of the topology of the non-zero one 
loop diagrams as reported in ( p,19|/]3~2l| / |3.22|) and Figs. |8| to [1C]. As noted in the previous 



section, the leading order contribution in each self energy, vacuum polarisation or vertex 
correction graph is obtained when only scalar gluons couple to quarks, and when all gluons 
are on shell, i.e. either soft or ultrasoft. Because for the Coulomb gauge, this eliminates 
any correction which is of the same order as the quark propagator or the Q^AqQ vertex 
(table |I|), the soft and potential quark self energy and all vertex corrections, Abelian or 
non- Abelian, are therefore zero to all orders in g and v : 

^Coulomb = x ^ ^Coulomb = j _ ^34) 

The coupling constant is hence renormalised only because the potential gluon vacuum 
polarisation in the Coulomb gauge is non-zero. A computation of the gluonic part which 
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makes again use of split dimensional regularisation can be found in eq. (41)]: 
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(3.35) 



The ghost part is zero because ghosts do not propagate. The contribution for light 
fermions is gauge independent and can hence be extracted from the Lorentz gauge result 
(|3.26 ). One arrives finally at 
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Not surprisingly, the NRQCD (3 function is extracted as (|3.32| ) using (|3.31| ). This par- 
ticular part of renormalising NRQCD in the Coulomb gauge may appear simpler than 
in the Lorentz gauges, but other parts suffer severely from technical problems like the 
non-transversality of the gluon propagator. 



4 Conclusions and Outlook 



After presenting the extension to NRQCD of a recently proposed explicit velocity power 
counting scheme for a non-relativistic toy field theory ]nj following Beneke and Smirnov's 
threshold expansion [IB], this article presented the computation of the NRQCD (3 function 
to lowest order in g and v in the Lorentz gauges and in the Coulomb gauge. It endorsed 
the relevance of a new quark representative and a new gluon representative in the soft 
scaling regime E ~ \p\ ~ Mv in which quarks are static and gluons on shell. HQET 
becomes a sub-set of NRQCD. The identification of three different regimes of scale for 
on-shell particles from the poles of the NRQCD propagators leads in a natural way to this 
regime, in addition to the well known potential one with on-shell quarks and instantaneous 
gluons mediating the Coulomb binding JR] and the ultrasoft one with bremsstrahlung glu- 
ons jy|]. Neither of the five fields in the three regimes should be thought of as "physical 
particles". Rather, they represent the "true" quark and gluon in the respective regimes 
as the infrared-relevant degrees of freedom. None of the regimes overlap. Using a rescal- 
ing technique |TJ|, |Tj], Sect. [] proposed an NRQCD Lagrangean which leads to the 
correct behaviour of scattering and production amplitudes p~5|. It establishes explicit 



velocity power counting, preserved to all orders in perturbation theory once dimensional 
regularisation is chosen to complete the theory. The reason is the non-commutativity 
of the expansion in small parameters with dimensionally regularised integrals and the 
homogeneity of dimensional regularisation. 

There is an intimate relation between the recently discovered threshold expansion |16] 



and this version of NRQCD which uses dimensional regularisation and has been developed 
over the past years [|14], [15], [18], [19[ . This was demonstrated at the outset of the compu- 
tation of the (3 function in Sect. |3.1| and has already led to fruitful exchange. As noted 
by Beneke and Smirnov threshold expansion provides an efficient way to derive the 
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NRQCD Lagrangean (|1.1|) from QCD, deepening our understanding of effective field theo- 
ries: The coefficients q, di, are obtained by computing QCD diagrams and sub-diagrams 
in the hard regime, q ~ M. Threshold expansion also shows how to systematise identify- 
ing the different kinematic regimes of NRQCD. Indeed, NRQCD was incomplete without 
the soft regime and its corresponding degrees of freedom as pointed out by Beneke and 



Smirnov fll6fl . One might therefore see threshold expansion and NRQCD as two sides of 
the same coin. The one expands loop integrals to establish a power counting for Feynman 
diagrams, while the other uses rescaling properties to classify the relative strengths of all 
vertices, allowing for intuitive interpretations of the derived vertex power counting rules 
on the level of the Lagrangean (see the end of Sect. |2.3| ). An effective field theory formu- 
lation can also be applied more easily to bound state problems. Finally, while threshold 
expansion has not yet been proven to be valid on a formal level, the existence of NRQCD 
as an effective field theory is guaranteed by the renormalisation group, e.g. [^, Chap. 8]. 

The calculation of the NRQCD (3 function in the Lorentz gauges to one loop order 
in Sect. ^| - albeit its result was easily anticipated - was non-trivial as it demonstrated 
that the power counting proposed is consistent also after renormalisation. Because of the 
splitting of quark and gluon fields in different representatives in the various regimes, this 
might be thought not to be straightforward. It also resolved a puzzle about the running 
of the coupling of gluons mediating the Coulomb interaction. The rescaling and power 
counting rules are gauge independent for a wide class of gauges including all standard 
gauges. Bare graphs and their corrections are of the same power in v. The soft regime is 
necessary not only to render the same f3 function as in QCD, but even to attain a result 
which is independent of the gauge parameter a at all and in which the renormalisation 
constants obey the Slavnov-Taylor identities. In the Coulomb gauge calculation of the (3 
function, the quark self energy and the vertex correction of the vertex coupling a quark 
minimally to a scalar gluon is easily shown to be zero to all orders in g and v. In the 
Coulomb gauge, the vacuum polarisation is therefore the only non-trivial renormalisation 
necessary for the (3 function. 

The derivation was greatly simplified by diagrammatic rules (Sect. |3.2| ) which follow 
from the homogeneity property of dimensional regular isat ion and allow one to systemat- 
ically recognise the majority of NRQCD graphs as zero to all orders in v just by drawing 
them, independently of the gauge and of details of the vertices involved. The rules apply 
equally well to threshold expansion and will prove fruitful in more rigorous investigations, 
also into NRQCD. Here, they were already used to prove scale conservation and the de- 
coupling of the ultrasoft, heavy quark in Sect. |3]3|. All topologies of one loop corrections 
to vertices involving one gluon and one quark were also classified, and the power counting 
of the leading order contributions determined. 

Dimensional regularisation in NRQCD has many advantages over other regularisation 
schemes in which scale-less integrals are non-zero. It preserves the symmetries and guar- 
antees the tree level power counting to be corrected at most by logarithms. It is well 
known that cut-off regularisation violates gauge symmetry and wants removal of power 
divergences in the renormalisation process in order to conserve the tree level power count- 
ing to all orders. When massless tadpoles do not vanish, the number of diagrams to be 
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computed in order to obtain even a simple vertex correction was seen to increase from 
one in dimensional regularisation to three or more. In general, it is not clear whether 
scale conservation holds in cut-off regularisation: The three different regimes require the 
introduction of at least three artificial cut-off scales to separate them and to cure ultra- 
violent divergences. Furthermore, the ultrasoft quark does not decouple, adding further 
diagrams. In the end, the answers of all well-defined regularisation schemes do agree, 
but dimensional regularisation provides an elegant way which furthermore often leads to 
analytic expressions. A disadvantage of dimensional regularisation is that in its standard 
formulation, it cannot be used to explore the non-perturbative sector of NRQCD. 

Finally, NRQCD with dimensional regularisation shows how to establish a power 
counting in any effective field theory with several low energy scales: The theory is first 
divided into one sector containing the particles relativistic at a given heavy scale and 
another sector containing the non-relativistic particles. In NRQCD, the former are light 
quarks, gluons and ghosts, the latter the heavy quark. The effect of all physics at larger 
scales is absorbed into the coefficients of the Lagrangean. Then, one identifies the combi- 
nations of scales in which particles become on shell by looking at the denominators of the 
various propagators. This gives the scaling regimes. Finally, the Lagrangean is rescaled 
to dimensionless fields in each regime to exhibit the vertex and loop power counting rules. 
A priori, all couplings obeying scale conservation are allowed. The diagrammatic rules 
developed in Sect. [O] still hold. Non-relativistic effective nuclear field theory is one ob- 
vious application with three low energy scales: the pion mass and production threshold, 
and the anomalously large scattering length of the two body system []. 

Among the points to be addressed further is an extension of the diagrammatic rules. 
An investigation of the influence of soft quarks and gluons on bound state calculations 



in NRQED and NRQCD is also important because - as seen at the end of Sect. |2T4 - 
their contribution at 0(g 4 ) and higher becomes stronger than retardation effects. To 
investigate the non-perturbative sector of NRQCD and to establish its equivalence to 
QCD in the continuum is another formidable task. 
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A Some Details on Split Dimensional Regularisation 



A.l Useful Split Dimensionally Regularised Integrals 

This appendix presents some useful formulae for non-covariant loop integrals using split 
dimensional regularisation as introduced by Leibbrandt and Williams |25j, see also the 
Appendix in [15]. In its results, split dimensional regularisation agrees with other meth- 
ods to compute loop integrals in non-covariant gauges, such as the non-principal value 
prescription p0[ , but two features make it especially attractive: It treats the temporal 
and spatial components of the loop integrations on an equal footing, and no recipes are 
necessary to deal with e.g. pinch singularities. Rather, it uses the fact that, like in ordi- 
nary integration, the axioms of dimensional regularisation [^, Chap. 4.1] allow to split 
the integration into two separate integrals: 



r d d k r 
J 72nV = J 



d a k d d - a k 



(2tt 



(27r) ff (2tt) 



d—a 



(a.i; 



Where applicable, split and ordinary dimensional regularisation of covariant integrals 
must hence agree once the limit a — ► 1 is taken with d still arbitrary. 
For the quark self energy, the integral 



d d q 



(2vr) a! (q 2 + ie) a q + a + ie 
is most easily computed by combining denominators as in HQET, 

d d q 



(A.2) 



2a / dA 



(2tt) £ 



q l + 2A(g + a) + ie 



>+i) 



(A.3) 



Now, the integrand is in a standard covariant form and can be computed using Minkowski 
space methods with the result 



1 



d d q 1 
(2n) d (q 2 + ie) a q + a + ie 



2 i T[| -a] F[2a + l-d] 
(4tt)5 r[a] 



(2a + ie 



.d-(2a+l) i7r(|-a) 



(A.l) 

By differentiating with respect to a, the following formula useful for the computation of 
the vertex corrections is derived: 



d d q 1 1 

(27r) d (q 2 + ie) a (q + a + ie) 



2 i T[f -a] T[2a + P-d] 
(4vr)i r[a] T{(3] 



(2a + ie) 



d-2a-/3 _i7r(f-a) 



(A.5) 

Another generalisation is easiest computed by first combining the denominators linear in 
go using Feynman parameters, 



1 



1 



q + a + ie q + b + ie 



dx 



1 



[q Q + xa + (1 — x)b + ie] 



(A.6) 
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and employing 
d d q 



1 



(27r) d {q 2 + ie) a q + a + ie q + b + ie 

-2 i r[f - a] T[2a + 1 - d] (2a + ie)^ 2 ^ 1 - (26 + ie)^- 1 
(4tt)^ r[a] 



(A.7) 



Note that in all three integrals ( A.4 / A"3 / A.7 ), a shift of g by an arbitrary value leaves 
the result unaffected. A scale in the momentum integration is only induced by the scales 
a, b of the energy integration. 



A. 2 Computation of Exemplary Vertex Corrections 

In order to illustrate the NRQCD power counting scheme further, this appendix outlines 
the computation of characteristic Abelian and non-Abelian vertex corrections in Sect. BTB 
The Feynman rules give for the Abelian correction to the QlA s Q Q s 



(t,p) 



(T',P') 



-i g R ) 3 ^ t b t a t c 



& d q —\8 



be 



1 - (1 -a 



4 



T + q T' + q 



(A.f 



(27r) d q 2 

(T-T',p-p') 

As usual, t b t a t b = (Cf — y)£ a - The integral containing q 1 in the numerator is simplified 
by re-writing it as 

rp rp/ rp rp/ 



So 2 



(T + q )(T + q ) 



(T + q )(T' + q ) (T + q ) (T> + q ) 



(A.9) 



and noting that the first integral is zero as it is without scale ( |3.10|) . The remaining 
integrals are using (|A.4j/|A.7|) 



^ R/ I C F t a re + finite, 

8tt 2 V 2 J 2 u 



(A.10) 



from which ( 3.28|) is read up. The Abelian correction to QiA p Q 



(T,p) 



(T-T',p-p') 



(T',P') 



. ., , , , / d d q —iS' 
-ig K f^t b t a t c ' 



be 



(2 



7T 



I - CI - a) q 4 



(A.ll) 



is inferred by replacing T — > T — and T' — > T' — g in (|A.10[) . Again, T 



2 / 

2M 



''Abelian 



IS 



extracted as quoted in ( p. 28 ). Note that - if one were to employ the equations of motion 
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for the outgoing particles - the integral would appear to vanish as no scale is present. 
This would mean that there is no Abelian vertex correction to the potential (or ultrasoft) 
scalar gluon vertex. In that case, the result for the NRQCD (3 function is not even gauge 
invariant. Since the diagram is the leading order contribution to the correction of the 
Q\>A P flQ p vertex (Fig. H), the application of the equations of motion is not legitimate. 



As noted in Sect. EM there are two Abelian corrections to Q\A S fiQ p . The leading one is 



(T,P) 



(T',p') 



i g R ) 3 ^ t b t a t c J 



d d q -i8 



be 



(2vr) d q 2 



:i-«)4 



Qo T + q 



. (A.12) 



The integration (|A.4j/SI7 ) gives the expected Abelian contribution fl3.28|) to the vertex 
normalisation (|3.27| ). 

The vertex correction already discussed in connection with the cutting rule, ( 3.13| ). 
is according to Sec. |3l] of order v 1 and hence not the leading contribution. Because 

Mv, one indeed finds 



T'~ |g ~ Mv 2 and T 



~ m ~ 



(T,P) 



(T',P') 



-i g R f ^ £ t b t a t c 



d^ — i5 



be 



(2vr) d q 2 



1-(1"«)4 



(T,p-p') - /2 

= ^ (Cf - y) *° ^ m + finite = P.C»! (A.13) 



The diagram will therefore only contribute to the renormalisation of the multipole expan- 
sion. Moreover, because it is sub-leading, one may use inside the integral the equations 

-*2 

of motion at leading order in the power counting, T = to set the outgoing potential 
quark on its mass-shell. The integral vanishes then, as no scale is present, see ( |3.13|) and 

(CT - 

The cutting rule suggests that the <3jAi,oQs correction should be related to the soft 
quark self energy: 



J : 

i g R ) 3 /i 3e t h t a f 



(T,P) ^ j -P- (T',p') 
(T-T'->0,p-p'-*0) 



(A. 14) 



d d q -i5 



be 



(2vr) rf q 2 



a] 



T + q T + q 



The dot at the cut-off vertex stands again as a reminder that the structure of the vertex 
will enter here, and one finds the vertex correction to be proportional to -^T, s (T,p) (|3.21 
Using ( |A.5| ), the vertex correction is again (|3.28| ). 
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Finally, two examples for the straightforward but tedious non-Abelian vertex corrections 
of Fig. ^] are given. First, one may calculate the non-Abelian correction to QlA u Q s 



(T,P) 



(T',P') 




g^(q) iG 0v (q) 



X 



d^ i i 

(2^p T + q ' q 2 q 2 

^Qo9^ + Qh9uq + Qu9^] ■ (A. 15) 



The Feynman gauge term, being independent of (1 — a), is zero as expected. For this 
diagram, even the terms proportional to (1 — a) 2 vanish, which can be shown to be finite 
in general and hence not to contribute to the vertex normalisation. In (|A.15| ), only the 
contribution leading in v is considered. Using the well known relation if abc t a t b = -| f 
for SU(iV), this can be collapsed down by re-writing numerators analogously to ( A. 9 ) to 



g 3 R ^Nt a (l-a)T J 



d d q 1 1 
(2^p T + q q~ 4 



T 2 



(A.16) 



One extracts the pole part of the non-Abelian vertex correction as ( 3.29Q from the result 
obtained using ( |A.4|) : 



i 9l A* 8 " N +a 3(1 - a) 



(An) 2 2 1 



T[e] + finite. 



(A.17) 



The leading piece of the non-Abelian correction to <5p^4 p ,oQp is extracted as 



(T,P) 



{T',P') 
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(T-T>,p-p>) 



■,2 n ,3s f a bCfbf c [ d d q i ig°»(q + AE) 

-ifla) 9k» f tt J—-—— (Q + AE) 2 



X 



(2n) d q q 2 (q + AEf 
x [2q g lu , + (AE-q) ft g vQ -(q + 2AE) v g lja ] , (A.18) 



where AE = (0,p' — p). The term proportional to (1 — a) contains the ultraviolet 
singularity and can be written as 



gi^Nt a (l~a) J 



1 



(2vr) d g 4 (q + AE)< 



q*-q*-2E-q) . 



(A.19) 



Standard formulae for dimensional regularisation yield as divergent part ( |3.29|) . The term 
proportional to (1 — a) 2 is ultraviolet finite and can be reduced to the integral 



& f | f (1 - af AE" J 



d d q ql 1 



{2n) d g 4 (q + AEY (4tt) 2 2 



i gl a 3£ N (1 — a) 2 

' — t a - - . (A.20) 
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